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Abstract 

In order to estimate roughness exponents of interface growth models, we propose 
the calculation of effective exponents from the roughness fluctuation a in the steady 
state. We compare the finite-size behavior of these exponents and the ones calculated 
from the average roughness {W2) for two models in the 2 + 1-dimensional Kardar- 
Parisi-Zhang (KPZ) class and for a model in the 1 + 1-dimensional Villain-Lai-Das 
Sarma (VLDS) class. The values obtained from a provide consistent asymptotic 
estimates, eventually with smaller finite-size corrections. For the VLDS (nonlinear 
molecular beam epitaxy) class, we obtain q = 0.93 ±0.01, improving previous esti- 
mates. We also apply this method to two versions of the ballistic deposition model 
in two-dimensional substrates, in order to clarify the controversy on its universality 
class raised by numerical results and a recent derivation of its continuous equation. 
Effective exponents calculated from a suggest that both versions are in the KPZ 
class. Additional support to this conclusion is obtained by a comparison of the full 
roughness distributions of those models and the distribution of other discrete KPZ 
models. 
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1 Introduction 



Surface and interface growtli processes attract much interest for tlieir appli- 
cations in solid state physics and material science [1,2,3]. They motivated the 
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study of many interface growth models, which led to significant advances in 
non-equilibrium Statistical Mechanics and related fields [2,3]. 

The simplest quantitative characteristic of a given interface is its roughness, 
also called the interface width, defined as the rms fluctuation of the height 
around its average position. The squared width, W2 = h? — h , is usually aver- 
aged over different configurations, and its scaling on time and length is used 
to characterize the growth process. For short times, in the so-called growth 
regime, the average roughness increases as (^2)3^ ^ ^"^^ i where /? is the growth 
exponent. For long times in a finite substrate, a steady state is attained, where 
the average roughness scales with the linear size L as 

{W2) - (1) 



with a called the roughness exponent (notice that steady state averages are 
denoted without further labels here, while indexes s or sat are frequently 
used in the literature to denote saturation). The connection of discrete or 
continuous growth models to certain growth equations is frequently based 
on accurate numerical calculation of those scaling exponents. However, this 
procedure may be quite difficult when huge finite-size or finite-time corrections 
are present in the scaling of {W2). 

Alternatively, one may investigate features of the full roughness distribution, 
P{w2), which is the probability density of the roughness of a given config- 
uration to lie in the range [-^2, ^2 -|- ^^2]. In recent years, the steady state 
distributions of several interface models have been analyzed [4,5,6,7] and they 
were shown to fit the scaling form 



Alternatively, the scaling relation 



can be adopted, where 



a ^ ^{W2^) - {W2f (4) 



is the rms deviation of the squared roughness. Despite the study of rough- 
ness distributions being useful for determining the universality classes of some 
growth processes, it is usually performed independently of the calculation of 



2 



scaling exponents (notice that the values of the exponents are not necessary 
to fit the scaled distributions). Indeed, finite-size corrections in the universal 
functions $ and \I/ are typically much smaller than those in the Family- Vicsek 
scaling relation [8], which gives Eq. (1) a limiting case. 

On the other hand, from Eqs. (3) and (4), we expect that the roughness fluc- 
tuation a in the steady state scales with exponent a as the average roughness: 

a ~ L'°. (5) 



Consequently, numerical estimates of a may also be used to estimate a, but, 
as far as we know, this was not done in previous works. The expectation of 
weaker finite-size corrections using a instead of (u;2) follows from the recent 
observation [7] that Eq. (3) provides a better data collapse of simulation data 
when compared to Eq. (2), for several growth models. 

In order to estimate a scaling exponent with accuracy, the standard method is 
to extrapolate effective exponents obtained from consecutive slopes of log-log 
plots (see, e. g.. Appendix A of Ref. [2]). Here, we will show that the effec- 
tive exponents obtained from a provide reliable estimates of the roughness 
exponent and, in some cases, varies with L slower than the estimates from the 
average width (w2)- This is illustrated with applications to three growth mod- 
els which belong to different universality classes, in one- and two-dimensional 
substrates. 

Subsequently, we will apply this method to calculate the roughness exponent 
of two versions of the ballistic depositon (BD) model [9] in two-dimensional 
substrates. Numerical works showed discrepancies in the estimates of that ex- 
ponent in d = 2 and a rigorous derivation of a continuous equation for BD [10] 
deviates from the Kardar-Parisi-Zhang (KPZ) equation [12]. However, when 
effective exponents of BD are calculated with cr, the asymptotic estimates of 
a are close to the best known KPZ values. The conclusion that this is the 
universality class of BD is reinforced by the comparison of the full roughness 
distributions. 

The rest of this work is organized as follows. In Sec. 2 we present the models 
which will be used to test the method for estimating roughness exponents and 
the results of these tests. In Sec. 3 we apply that method to two versions of BD 
in (i = 2 and analyze their roughness distributions. In Sec. 4 we summarize 
our results and conclusions. 
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2 Tests of the method in models of the KPZ and the VLDS classes 



The first model analyzed here is the RSOS model [13], in which the incident 
particle can stick at the top of a column only if the differences of heights of 
all pairs of neighboring columns do not exceed AHmax = 1 after aggrega- 
tion. Otherwise, the aggregation attempt is rejected. The second model was 
proposed for etching of a crystalline solid by Mello et al [14] : at each growth 
attempt, a randomly chosen column i, with current height h{i) = ho, has its 
height increased of one unit {h{i) — > ho + 1), and all neighboring cohimn whose 
heights are smaller than ho grows to ho (this may be called the growth version 
of the etching model [11]). 

In the limit of large lengths and long times, these models are represented by 
the KPZ equation [12] 

^^U2V'h + X2{Vhf + rj{x,t). (6) 



Here, 1^2 and A2 are constants and 77 is a Gaussian noise with zero mean and 
variance [r] {x, t) r] (x', X — X ''^5{t-t'), where d is the dimension 

of the substrate. In = 1, the exact roughness exponent is a = 1/2 [12,2], and 
md = 2 the best numerical estimates arc around 0.38 — 0.39 [15,11]. Numerical 
estimates of (W2) for the RSOS and etching models \nd—l provide estimates 
of exponents which converge to the exact values with negligible finite-size 
corrections. However, that is not the case in = 2, particularly for the etching 
model [11], consequently our tests will focus this case. 

The widths (u'2) and their fluctuations a were calculated at the steady states in 
two-dimensional substrates up to L = 256 for the RSOS model and up to L = 
512 for the etching model. Finite-size estimates of the roughness exponents 
are given by 

l ln[(^2) {L)/{w2){L/2)] 



and 

Un[a{L)/a{L/2)] 

"^^^^^ = 2 hr2 ■ 

The presence of significant corrections to Eqs. (1) and (5) is reflected in a 
size-dependence of (L) and (L). These effective exponents are shown in 
Fig. 1 for the RSOS and the etching models as a function of 1/L and 2/L, 
respectively. 
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An asymptotic estimate near 0.39 is consistent with the trend of the RSOS 
data as L oo. For this model, small differences in the estimates obtained 
from (^2) and a are observed, as well as weak finite-size effects. The linear fit 
of a^, (L) shown in Fig. 1 suggests a ~ 0.39. 

However, the finite-size dependence of the data for the etching model is re- 
markable: the effective exponents vary between 0.34 and 0.38 in Fig. 1. For 
the largest lengths L, (L) approaches the asymptotic region a ~ 0.39 faster 
than aw{L), as illustrated by the curves through the data points in Fig. 1. 
It suggests that, for large L, the corrections to the scaling of a are smaller 
than the corrections to (^2)- Although the fits of the data in Fig. 1 do not 
represent systematic extrapolation procedures, our interpretation is reinforced 
by the fact that the extrapolated (L) for the etching model is very close to 
the value predicted by the RSOS model. 
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Fig. 1. Effective exponents (L) (empty symbols) and (L) (full symbols), for 
the RSOS model (squares) and the etching model (triangles) in two-dimensional 
substrates, as a function of 1/L or 2/L. [L) data are slightly shifted to the right 
to avoid superposition of data points. Error bars are compared to the size of the 
points, except for the largest lengths, in which the uncertainties are around 0.005 
for both models. The lines connecting the data points (except (L) data for the 
RSOS model) were drawn to guide the eye. 

For the etching model, we were not able to perform reliable extrapolations of 
ayj (L) or fto- (L) as functions of l/L^, with any exponent A > (the value 
A = 1 is used in Fig. 1). This is related to the large deviations from the 
asymptotic scaling in small lattices, typically L < 50. It contrasts to Ref. [11], 
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where a extrapolation of a^, (L) up to L = 1024 was possible with A ~ 0.5. 
However, here we are working with data for L < 256 because there are very 
large error bars in the estimates of a for larger lattices. With a small number 
of data points, even more sophisticated extrapolation methods [16] are not 
able to provide better estimates than previous work. 

An extended version of the conserved RSOS (CRSOS) model is also suitable 
for this test. In this model, if the aggregation at the column of incidence dis- 
obeys the condition AH max = 1) then the incident particle executes a random 
walk among neighboring columns until finding a column in which it can aggre- 
gate [17]. The CRSOS model belongs to the class of the nonlinear molecular- 
beam epitaxy equation, or Villain-Lai-Das Sarma (VLDS) equation [18,19] 



^ = j,,W'h + X,V'iVhf + rjiS, t), (9) 



where U4 and A4 are constants. Our extended version of the CRSOS model is 
different from the original one [20] but has the same symmetries and, conse- 
quently, belong to the same universality class. 

In Fig. 2 we show the effective exponents of this model obtained from simula- 
tions up to L = 256, with linear fits of both data sets. The evolution of (L) 
and aw {L) suggest a ~ 0.93 asymptotically, but the difference between (L) 
and the asymptotic value is clearly much smaller. This is consistent with the 
different slopes of the hnear fits shown in Fig. 2 (smaller slope for fitting 
a. (L)). 

This application provides additional support to the conclusion that a < 1 in 
the VLDS class in li = 1 [21] (a = 1 is obtained from one-loop renormaliza- 
tion). The estimate a = 0.93 ± 0.01 follows from the extrapolation of those 
data. It improves our previous estimate a = 0.94 ± 0.02, obtained in Ref. 
[17], where effective exponents were calculated from the second and fourth 
moments of the height distribution in sizes up to L = 1024. 

The above analysis, performed with models in different universality classes and 
with noticeable (but not huge) finite-size corrections to Eq. (1), shows that 
roughness exponents may be estimated from a with comparable or higher ac- 
curacy than those obtained from W2, suggesting its application to controversial 
situations. 
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Fig. 2. Effective exponents (L) (empty squares) and a^j (L) (full squares), for the 
CRSOS model in one-dimensional substrates, as a function of l/L. Error bars are 
of the size of the data points. The dotted and dashed lines are linear fits of each 
data set. 

3 Application to ballistic deposition 



In the simplest version of the ballistic deposition model, particles are released 
from a randomly chosen position above a (i-dimensional substrate, follow tra- 
jectories perpendicular to the surface and stick upon first contact with a near- 
est neighbor occupied site [8,9]. The resulting aggregate is porous and has a 
rough surface. Several applications of the BD model or its extensions to real 
growth processes were already proposed, which justifies the present analysis 
(see e. g. recent applications in Refs. [22,23]). 

The mechanism of lateral aggregation suggests that BD is in the KPZ class. 
However, several numerical works on this model showed discrepancies between 
the estimated exponents and the KPZ values [8,24,25,26]. In a recent work, we 
have shown that the KPZ values are obtained with good accuracy in d = 1 if 
effective exponents (L) are properly extrapolated [27]. In d = 2, the effec- 
tive roughness exponents rapidly vary with L, but they are still far from the 
region [0.38,0.40] up to L = 512 [27]. An apparent consistency with the KPZ 
value is obtained in ci = 2 only after assuming of a constant correction term in 
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Eq. (1) (the intrinsic width), but error bars arc very large [11]. Consequently, 
from numerical work, there is no clear evidence that BD is in the KPZ class 
md = 2. 



Recently, Katzav and Schwartz derived a continuous equation for another ver- 
sion of BD which differs from the KPZ equation [10] . They analyzed the BD 
model with next-nearest neighbor aggregation, hereafter called BDNNN, in 
which the incident particle sticks upon first contact with a nearest neighbor 
or with a next-nearest neighbor occupied site. In d — 1, they showed that the 
squared gradient in Eq. (6) is replaced by the absolute value of the gradient 
in the continuous equation corresponding to BDNNN. This is suggested as 
a possible reason for the slow convergence of the scaling exponents to their 
asymptotic values. In o? > 2 the situation is more complicated: the only contri- 



ong which the 
. This ac- 



dh 



dXr, 



bution to the local growth rate comes from the direction a 
height gradient is maximum, and that contribution is + 
counts for the maximal growth ingredient of the model. However, the presence 
of the lattice axis in that equation contrasts with the rotation symmetry of 
the KPZ equation. Thus, also from analytical work, there is no a priori reason 
to expect BD to be in the KPZ class [10]. 

This scenario motivates the extension of our method to estimate roughness 
exponents of BD and BDNNN in d = 2. Simulations of both models were 
performed in lattices up to L = 512, where W2 and a can be obtained with 
reasonable accuracy. The number of different realizations used to estimate 
those quantities were near 10^ for L < 128, 10''' for L = 256 and 2 x 10^ for 
L = 512. 



In Fig. 3a wc show (L) and «o- {L) for the original BD model. (L) rapidly 
varies with L in the whole range of our data and attains a value near 0.28 in 
the largest lattice, which is much smaller than the estimates for the KPZ 
class. This complex finite-size behavior was already pointed out in Ref. [27] 
and, consequently, no reliable extrapolation of a could be extracted from those 
data. On the other hand, a^j (L) is in the range [0.36, 0.37] for the largest L. 
As 1/L — > 0, Fig. 3a indicates that a„ (L) converges to an asymptotic value 
between 0.36 and 0.4. 



In Fig. 3b we show the effective exponents for the BDNNN model. We notice 
that (L) is even farther from the KPZ estimate in this model. However, 
the trend of cto- (L) as 1/L — > suggests an asymptotic value between 0.35 
and 0.4, although finite-size effects seem to be much stronger than those of 
the original BD model. 

Extrapolations of a^, (L) and (L) for BD and BDNNN as functions of 1/ L^, 
varying the exponent A along the same lines of Ref. [11], were also tested, but 
no significant improvement of the above mentioned estimates was obtained. 
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Fig. 3. Effective exponents (L) (empty symbols) and (L) (full symbols) for: 
(a) the original BD model in d = 2; (b) the BDNNN model in d = 2. The lines 
connecting the data points were drawn to guide the eye. 

However, those estimates are consistent with the best currently avaiable values 
for the KPZ class, which lie in the range [0.375,0.396] [15,11], which strongly 
suggests that BD and BDNNN are also in that class in d = 2. 

Another important test is a comparison of the full roughness distributions of 
these models and that of other KPZ models. The distribution for the RSOS 
model, which shows negligible finite-size effects, can be used as a representative 
of the KPZ class [7]. In Fig. 4 we show this distribution, scaled according to 
Eq. (3), and the scaled distributions for the BD and the BDNNN models in 
L = 256 (data for L = 512 was not used here due to their lower accuracy, 
mainly in the tails of the distribution). The collapse of all data into a single 
universal curve reinforces our conclusion that both models are in the KPZ 
class in d = 2. 
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Fig. 4. Scaled roughness distribution at the steady states of the RSOS model 
(solid curve), the BD model (squares) and the BDNNN model (triangles), in 
2 + 1-dimensional lattices of length L = 256. 

A quantitative support to the universality of the scaled distributions of Fig. 4 
is given by the estimates of the skewness and kurtosis of those curves. For the 
RSOS and the etching models, S = 1.70±0.02 and Q = 5.4±0.3 were obtained 
in Ref. [7]. Here, for BD we obtained S = 1.71 ± 0.02 and Q = 5.4 ± 0.1, and 
for BDNNN we obtained S = 1.67 ± 0.03 and Q = 5.2 ± 0.3, both in good 
agreement with the values of the other KPZ models. 



4 Conclusion 



We proposed the calculation of roughness exponents from the scaling of rough- 
ness fluctuations a in the steady states. Tests for two models in the 2 + 1- 
dimensional KPZ class and for a model in the 1 + 1-dimensional VLDS class 
were performed and showed that effective exponents obtained from that quan- 
tity provide comparable or better asymptotic estimates than the average width 
(^2)- We also applied this method to two versions of the ballistic deposition 
model, the original one (BD) and the next-nearest neighbor one (BDNNN), 
in order to clarify a controversy on its universality class raised by numerical 
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results and an exact derivation of its corresponding contimioTis equation. Ef- 
fective exponents calculated from a suggest that the asymptotic class is KPZ, 
while those calculated from (W2) show much more complex finite-size behav- 
ior. Additional support to this conclusion is obtained from the collapse of the 
full roughness distributions of those models into the distribution of the RSOS 
model. 

At this point it is important to recall that higher moments of the height 



distribution, such as the fourth moment (W4) = i^h — h y, are also frequently 

used to estimate the roughness exponent. Their accuracy is usually lower than 
that of (^2) and, consequently, do not improve the estimates of a (see e. g. 
Ref. [11]). This quantity, however, must not be confused with cr^, which is the 
fluctuation of the global roughness in the steady state. 

The successful application of the above method to estimate the roughness ex- 
ponent, possibly combined with the analysis of roughness distributions, sug- 
gests its extension to other models in which remarkable finite-size corrections 
in the effective exponents are observed. 
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